Point Contact Spectroscopy of Superconducting Gap Anisotropy in Nickel 

Borocarbide Compound LuNi 2 B 2 C 
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Point contacts are used to investigate the anisotropy of the superconducting energy gap in 
LUNI2B2C in the ab plane and along the c axis. It is shown that the experimental curves should 
be described assuming that the superconducting gap is non-uniformly distributed over the Fermi 
surface. The largest and the smallest gaps have been estimated by two-gap fitting models. It is 
found that the largest contribution to the point-contact conductivity in the c direction is made by 
a smaller gap and, in the ab plane by a larger gap. The deviation from the one-gap BCS model is 
pronounced in the temperature dependence of the gap in both directions. The temperature range, 
where the deviation occurs, is for the c direction approximately 1.5 times more than in the ab 
plane. The Y parameter, allowing quantitatively estimate the gap anisotropy by one-gap fitting, in 
c direction is also about 1.5 times greater than in the ab plane. Since it is impossible to describe 
satisfactorily such gap distribution either by the one- or two-gap models, a continuous, dual-maxima 
model of gap distribution over the Fermi surface should be used to describe superconductivity in 
this material. 

PACS numbers: 63.20.Kr, 72.10.Di, 73.40.Jn 
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INTRODUCTION 



The goal of this study was to investigate the anisotropy 
of the energy gap in nickel borocarbide superconductors 
i?eNi2B2C. The crystallographic structure of these com- 
pounds resembles to some extent the structure of HTSC 
materials. 1 It consists of alternating Re-C and M2-B2 
layers (Re is a rare earth metal). There is a point of view 
that the electron properties of i?eNi2B2C compounds in 
the normal state are isotropic because of the strong car- 
bon bond along the tetragonal c axis^ This is supported 
by the temperature dependence of isotropical resistivity 
in the single crystal of YNi2B2C£ On the other hand, 
there is also data about a substantial anisotropy in these 
compounds^ In the superconducting state the experi- 
mental data are also contradictory. The upper criti- 
cal fields for H parallel and perpendicular to the c axis 
and the derived superconducting parameters do not show 
any anisotropy for the YNi2B2C single-crystal samples 
in agreement with magnetization and torque magnetom- 
etry measurements, but a small anisotropy is observed 
for the LUN12B2C single crystals*^ Meanwhile, a very 
strong anisotropy of the superconducting energy gap was 
reported for LUNI2B2C in the ab plane. This conclusion 
was based on the observation of delocalized quasiparticles 
in thermal conductivity at very low temperatures^ Ac- 
cording to this analysis, the smallest gap differs at least 
ten times from the gap in other directions. Recently, evi- 
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dence for the presence of nodes along (100) direction was 
provided by the field-angle thermal conductivity, 7 field- 
angle heat capacity 8 and ultrasonic attenuation measure- 
ments of YNi 2 B 2 CiS 

In accordance with the experimental data an 
anisotropic superconducting gap function was proposed 
in Ref.[lfJ (see, also, Fig. 52 from the recent reviewii). In 
that model the gap has nodes along (100) directions and 
attains the maximum values along (110) directions. In 
the case of s + g symmetry, the elastic scattering leads 
both to a decrease of T c and to an "isotropization" of the 
gap with vanishing nodes. 

We do not know any direct measurement of gap 
anisotropy involving the Andreev reflectionii 2 . The STM 
tunnelling measurements at 4K report a gap of 2.2 meV 
along the c axis, yielding a too low ratio 2A/fcT c = 3.2 
for LuNi2B2C4ii4 By contrast, there are point-contact 
measurements of the gap in the ab plane (e.gi^^), which 
yield 2A/kT c = 3.7-^3.8. The point-contact spectroscopy 
method in the latter measurements cannot, however, pro- 
vide angular resolution of the gap anisotropy much better 
than 7r/2. 

In the present investigation we have found that both 
in the ab plane and in c direction the experimental 
point-contact spectra cannot be fitted satisfactory with 
the one-gap theoretical curve, even when broadened by 
an adjustable broadening parameter. Fitting the ex- 
perimental spectra both at a small bias (eV < A) 
and at a bias larger than A forces us to use at least 
the two-gap fitting curve like the one for the two-band 
superconductori 12 ! 17 But even the two-gap fitting with 
the two proper broadening parameters was not good 
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enough to approximate the experimental characteristic 
at the middle energy region (eV ~ A). Moreover, for the 
two-band model one could expect transitions of Cooper 
pairs between the bands to introduce an additional de- 
pairing factor, which is analogous to other depairing fac- 
tors like magnetic impurities, strong magnetic fields, etc. 
In such case the superconducting order parameter and 
the energy gap should differ from each other, which could 
be accounted for by using Ref.^l Unfortunately even 
using the latter theory satisfactory fits could not be ob- 
tained. Therefore we suggest that only a continuous gap 
distribution with the two energy gap maxima around the 
Fermi surface could satisfactorily approximate the exper- 
imental spectra at low temperature (T << T c ). We have 
determined tentatively the value of these maxima and 
the range of gap distribution which corresponds well to 
the recent STM measurements^ 



II. EXPERIMENTAL TECHNIQUE 

The point contact measurement was performed on sin- 
gle crystal LuN^E^C with T c ~ 16.9 K grown by Can- 
field and Bud'ko using a flux method^ Geometrically, 
our crystal was a thin (0.1 ~ 0.2 mm) plate with the c 
axis perpendicular to its plane. The single crystal sur- 
face always contains quite a thick layer in which super- 
conductivity is either absent or strongly suppressed. To 
perform measurement in the ab plane, the crystal is usu- 
ally cleaved and the point contact is made between a 
metallic counterelectrode and the cleaved surface. It was 
technically problematic to produce a cleavage perpendic- 
ular to the c direction. In this case the crystal surface 
was cleaned with a 10% HNO3 solution in ethanol. As 
the measurement in the ab plane shows, both cleavage 
and etching yield identical results. 

The point contact in the ab plane is fabricated between 
an edge of the silver electrode and a freshly cleaved (or 
etched) corresponding facet of the single crystaliSi The 
deviation from the perpendicular to the c direction might 
amount about 5 -7- 10°. We do not know a priori along 
which of the in-plane directions the contact is obtained. 
But since we used a selection rule to choose the highest 
observable superconducting energy gap with the largest 
nonlinearity at the gap double-minimum structure, the 
contact axis is presumably along (110) directions, where 
the in-plane gap is maximum. To produce contacts in 
the c direction the traditional "needle-anvil" geometry 
is used. The radius of the needle is about 1-3 micron. 
The temperature was measured using a special cryogenic 
insert (its close analogue is described elsewhere^) 

The point contact resistances varied typically from sev- 
eral Ohms to tens of Ohms. For detailed investigations 
we chose the point contacts with the highest permissible 
tunneling which was controlled against the differential 
resistance maximum at zero bias and by the maximum 
nonlinearity in the dV/dl-double- minimum region, cor- 
responding to the intact superconducting surface under 



the contact. Unfortunately, a complete set of curves in 
the whole T m i n — T c interval was obtained only for a few 
contacts. Because of their rather high resistance and the 
long-duration (over 10-12 hours) of the measurements, 
many of the contacts were broken down. The tempera- 
ture interval ranges from the minimal available T (~1.5 
K) to T,. 



III. RESULTS AND DISCUSSION 

Some curves of the first derivatives dV/dl measured at 
different temperatures on LuNi2B2C-Ag point contacts 
along c axis and in the ab plane are shown in Figs. 11121 



48 

46 

44 

42 

§ 40 

5 38 
> 

-° 36 

34 
32 



1 1 1 1 1 




\ /^T, K - 


; \ 


/ 1.47 - 




\ 1 513 : 


;LuNi 2 B 2 C-AgA ' 


1 80 ' 


"c-direction !u 


J' 12.0 ~ 




J 15.0 ; 




3 175 _ 


.... 1 1 


i i . . . . 



-20 -10 10 
Voltage (mV) 



20 



FIG. 1: Differential resistance of LuNi2B2C-Ag point contact 
at different temperatures. The contact axis is parallel to c. 
Not to overload the figure, we give only several representative 
curves. 

The measured curves were symmetrized dV/dI sym = 
1/2 [dV/dI(V)+dV/dI(-V)] and normalized to the nor- 
mal state at T > T c , except in Figs 1,2, where the raw 
data are shown. Each curve contains several hundred of 
experimental dots and at the scale shown in presented fig- 
ures, all the curves have negligible noise. In other figures 
experimental curves are shown by dots whose number is 
decreased in order to be discerned. Statistically (several 
tens of contacts were examined for each direction), the 
distance between the minima in the first derivative dV/dl 
characterizing the average value of the gap was 12-15% 
larger in the ab plane compared to that in the c direction. 
This is a very crude estimation of the c <-> ab anisotropy. 

The theoretical predictions and the experimental re- 
sults were compared using two approaches. First, a 
model was applied, which describes electrical conductiv- 
ity of pure S-c-N point contacts in the presence of an 
arbitrarily transparent potential barrier at the boundary 
between the metals. This model allows for the finite life- 
time of the Cooper pairs.™ The I — V characteristics of 
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FIG. 2: Differential resistance of LuNi2B2C-Ag point contact 
at different temperatures. The contact axis is along the ab 
plane. See the last sentence in the previous caption. 



the point contact are described as follows: 
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is the 



transparency of the potential barrier, and Z is the tun- 
neling parameter. The value u = e/ A e (A e - is a complex 
energy-independent order parameter of the superconduc- 
tor), which can be found from the equation^ 



where A is the BCS order parameter obtained from 
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and 7 = 1/rgA is the pair breaking parameter (rs is 
the mean free time during, e.g., spin-flip scattering at 
impurities). When magnetic impurities are absent, t$ 
tends to infinity and Eq. (1) coincides with the equation 
for current in Ref.|3 The energy gap Ao and the order 
parameter A are related as follows: 



A„ = A(l-7 2 / 3 



3/2 



(4) 



The other approach was based on the generalized 
Blonder-Tinkham-Klapwijk (BTK) models commonly 
used to describe S-c-N point contacts. The model al- 
lows for the finite lifetime of quasiparticlcs t—Tl/T deter- 
mined by inelastic scattering, which leads to the broad- 
ening of the density of states in the superconductor. For- 
mally, according to the theorjii& the BTK-based results 
are obtained under the condition of strong pair breaking 
(|u|^>l). Therefore, in the strict sense, the generalized 
BTK model contains no gap. For any infinitesimal broad- 
ening parameter V, at T — the density of states near 
the Fermi surface is non-zero. In theory^ the order pa- 
rameter is a quantity analogous to the pseudogap in the 
generalized BTK model. In the following we use terms 
" order parameter" and 7 for the model of Refll8l while 
the "gap" and T for the commonly used BTK model^ 
except special cases, where we refer to Ao related to the 
models 

Now we consider more closely the methods of fitting 
the theoretical curves to the experimental results. The 
iteration method commonly used in this case is quite 
good for small broadening. However, when the broaden- 
ing T is comparable with the gap (or 7 > 0.3), the results 
thus obtained are ambiguous and often dependent on the 
starting A, T(or 7) and Z values. In our calculation, we 
therefore used the technique of coordinate descent with 
a postponed solution. 

First, we specified an interval in which A is searched at 
a given temperature. The interval was then subdivided 
into equidistant parts: Ai, A2, ...A n . 7 and Z were fit- 
ted for each Aj . The procedure was as follows: after each 
step of 7 -fitting, a complete Z-fitting was performed and 
then the next step of 7-fitting was considered. Before 
each calculation of the average r.m.s. deviation F(A^), 
the amplitudes of the fitting and the experimental curves 
were made equal by multiplying the y-coordinate of the 
fitting curve by a scale-factor S and then shifting it along 
the y-axis by an amount B. The values of B and S were 
chosen to minimize the deviation F(Aj). The standard 
algorithm for determination of B and S is considered, for 
instance, in Ref.|2^. As a result, for each Ai we found 7$ 
and Zi at which the difference between the shapes of the 
theoretical and experimental curves characterized by the 
r.m.s. deviation F(Aj) was the smallest one. The same 
holds for T in the BTK fitting models The calculation 
for some temperatures is shown in Fig. [3] It is evident 
that at T c the technique cannot ensure unambiguous re- 
sults. At 14.5 K, the curve F(Aj) is practically horizon- 
tal. The determination of the A-value should take into 
account other factors as well. For example, it is impor- 
tant that the tunneling parameter Z and the scale factor 
S should be invariant. The dependences S(A) and F(A) 
at T = 7.5 K are shown in Fig.0] It is seen that the 
minimum F corresponds to a certain S-value. The BTK 
calculation of the corresponding dependences at the same 
temperature yields broader F(A) curves. Nevertheless, 
the minima positions in these curves are quite definite. 

At low temperatures the shapes of the theoretical and 
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FIG. 3: Averaged r.m.s. deviation F characterizing the ex- 
tent of shape discrepancy between the theoretical and exper- 
imental curve as a function of A. The minima in the curves 
correspond to the best agreement of theory and experiment. 
The data shown are for LuNi2B2C-Ag point contacts in the 
c direction at different temperatures. Calculation based on 
Ref.E 
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FIG. 5: Temperature dependence of the shape discrepancy 
between theoretical and experimental curves calculated ac- 
cording to Ref.^^ and within the generalized BTK model— 
for LuNi2B2C-Ag point contacts in the c direction. 
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FIG. 4: A-dependences of the averaged r.m.s. deviation F 
(see Fig-EJ an d the scale factor S at T — 7.5 K. Calculation 
by equations from Ref.llSt S is a factor used to divide the 
j/-coordinate of the theoretical curve to match its amplitude 
with that of the experimental curve. 



FIG. 6: Approximation of the experimental curve with one- 
and two-gap models for the LuNi2E>2C-Ag point contact in 
the ab plane. The dotted straight lines marked with an arrow 
indicate the kink in the experimental curve which is exagger- 
ated by the 2-gap theoretical fit. 



the experimental curves were in rather poor agreement 
for all of the borocarbides investigated (i?e=Er, Dy, Tm, 
Lu) irrespective of the point contact orientation. The 
temperature dependence of the r.m.s. deviation in shape 
in the c direction is shown in Fig.[S]for both models>i2i2£ 
It is seen that the best agreement is achieved above 7 K. 
At T <7 K, the largest discrepancy between the shapes is 
observed in the region of extrema. The central maximum 
of the differential resistance in the experimental curve is 
considerably narrower than that in the best fit;2£ whereas 
the minima are shifted to higher bias (Fig.[SJ). 

Some experimental curves have a kink marked with an 



arrow in Fig.0 This shape is typical of a superconduc- 
tor with two gaps of close energies. The temperature 
dependence of the broadening (pair breaking) parame- 
ter also supports the existence of two gaps. Since Lu is 
a nonmagnetic material, it is natural to expect a negli- 
gible broadening (pair breaking) parameter for the Lu- 
based borocarbide point contacts. However, that is not 
the case, as is evident from the temperature dependences 
(Fig.JTJ of the broadening (pair breaking) parameters ob- 
tained within the two models (Ref .0 and BTK2£) . 

Note that 7 and T are different quantities. In the 
generalized BTK model 2 - 5 V is independent of the gap 
and has the gap dimensions; whereas, in the model of 
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FIG. 7: Temperature dependences of the broadening and pair 
breaking parameters for the LuNi2B2C-Ag point contacts in 
the ab plane (a), and along the c axis (b). The index "max" 
marks the "wings" fitting (see text). "Wings" fit for ab- plane 
gives 7=0. 



Ref.lTsl 7 is described as 7 = 1/rgA, and, for 7 < 1, 7 
is related to the gap Ao and the order parameter A as 
7 = [1 — (A0/A) 2 / 3 ] 3 / 2 . As the temperature is increased, 
both r and 7 decrease, and, at temperatures >7.5 K 
for the ab-plane and >11 K for the c-direction, tend to 
zero. Meanwhile, the one-gap fit starts to approximate 
the experimental curve more and more accurately. In this 
case the one gap fit with the broadening (pair breaking) 
parameter is a tool of describing a certain average gap 
(order parameter). 

The largest and smallest gaps can be estimated by fit- 
ting different portions of the experimental curve. In this 
case, to match the amplitudes of the fitting and the ex- 
perimental curves, we used another method. Namely, the 
y-coordinates at the central maximum and near the min- 
ima are scaled to coincide (Fig.[SJ). The best agreement 
in the central maximum region can give us an estimate 
of the smallest gap (curve 1), while the best agreement 
in the "wing" region yields the largest gap (curve 2). By 
"wings" we imply the portion of the experimental curve 
at biases higher than the energy gap double-minimum 
structure in the dV/dl characteristic, namely at eV > 3. 
The estimates however are rather rough. 

Further, we follow the method of approximation of I — 
V characteristics, which was calculated for the S — I — N 
contacts of MgB2 within a two-band model in Ref.llTi 



There, the total conductivity of the tunnel contact is 
a sum of the n— and a— band conductivities analyzed 
by applying the BTK model. To describe the resulting 
curve, we also used a model of two independent parallel- 
connected point contacts with different gaps whose con- 
ductivities are additive. The contributions of these con- 
ductivities account for the part of the Fermi surface con- 
taining a particular gap. Thus, for the two-gap model an 
experimental curve is fitted by the following expression: 



dV 
dl 



§(A 1 , 11 ,Z)K+§(A 2 , 72 ,Z)(1-K) 



(•5) 



with a proper choice of the coefficient B. Here, the co- 
efficient K reflects the contribution of the part of the 
Fermi surface having the gap Ai, S is the scaling fac- 
tor discussed for the one-gap approximation. To obtain 
the best agreement with the experiment, the parameters 
used in this expression are allowed to differ from those 
found for fitting of the separate portions of the exper- 
imental curve. It appears that when the contribution 
of the smaller gap prevails, the "wings" fit gives smaller 
values. In this case, the best result is achieved with the 
two-gap fit (Fig. 8, curve 4, asterisks) using the param- 
eters for curves 1 and 3, the latter is marked with dots. 
Correspondingly, if the larger gap prevails, the central- 
maximum fit gives higher values. In Fig. 6 an example of 
the two-gap fit (asterisks) for the aft-plane is shown. 

Although the two-gap fit shows much better agree- 
ment in the regions both of the central maximum and 
the "wings" (Figs. 6, 8 asterisks), it cannot provide a 
complete description of the experimental curves, espe- 
cially at the double-minima structure of dV/dl, and in 
the c direction (Fig.|SJ). Most likely, this is because a 
continuous distribution of the gaps. 

The previous estimation of point-contact parameters in 
the ab plane and in the c direction was made within one- 
gap and two-gap models using the equations of Ref.ll8l 
and BTK 25 . The results are presented in the Table. Ai, 
A2 are the highest and lowest order parameters (the same 
as BTK gaps^ 5 -) from the two-gap fit; A01 and A02 are 
the energy gaps corresponding to Ai, A2; 71, 72 are the 
pair breaking parameters of the modeli^ for the largest 
and smallest gaps. (The same holds for the widening 
parameters Ti, T 2 in the BTK model. 25 ) Independently 
determined Z's arc the tunneling parameters, which are 
the same for the largest and smallest gaps and are thus 
quite selfconsistent. In the Table the order parameters 
(gaps) based on the two-gap fit appear along with the 
contributions to the total conductivity from the Fermi 
surface region with the corresponding gap. For example, 
for the c direction contact, the contribution to the con- 
ductivity is 40% from the region with an order parameter 
of 3.2 mV and 60% from the region with 1.82 mV. 

It is seen that within both the modelsjiS*^ 5 . the dom- 
inant contribution to the total conductivity in the c di- 
rection is made by the Fermi surface region where the 
order parameter is lower. The higher order parameter 
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FIG. 8: The estimates of the smallest and largest gaps for the 
LuNi2B2C-Ag point contact in the c direction at T = 1.47 K. 
Curve 1: A = 1.82 mV (A = 1.65 mV), 7 = 0.016, Z = 
0.558; curve 2: A = 2.97 mV (A = 2.693 mV), 7 = 0.016, 
Z = 0.553; curve 3: A = 3.2 mV (A = 2.9 mV), 7 = 
0.016, Z — 0.558. The best two-gap approximation (curve 
4, asterisks) is achieved with the parameters of curves 1 and 
3. The relative contribution to the total conductivity is 60% 
from the smaller gap and 40% from the larger gap (see the 
Table). The experimental curve is shown with solid dots. 



makes the dominant contribution in the ab plane. On 
the average this correlates with gaps in ab and c direc- 
tions predicted by the (s + g) model— 

Although the part-by-part fitting is rather rough, we 
tried this procedure for tracing the temperature depen- 
dences of the smallest and largest gaps. However, more or 
less definite values were obtained only at the lowest tem- 
peratures because the smallest gap is estimated within a 
relatively small part of the experimental curve (near zero 
bias). As for the largest gap, its temperature dependence 
was traced up to the moment when the calculations over 
the entire curve and over the "wings" start to give similar 
results. 

The temperature dependences of the scale factor S in 
the c direction are shown in Fig.|5| They were found 
for the largest gap ( "wing" fitting is marked with "max" 
throughout) and for the average gap (entire-curve fitting) 
by the equations of Ref.0 and BTK— The tempera- 
ture dependences of the tunneling parameter Z are shown 
for both models. To avoid crowding in the figure, only 
the Z-values obtained from the entire-curve fitting are 
shown (the "wing" - fitting results are much the same). 
Although the fitting was focused on the coincidence of 
theoretical and experimental dV/dI(V)- curves and no 
parameters were specially fixed, we can state that Z and 
S vary very little in the whole temperature interval. At 
T = 9 K the "wing" and entire-curve fits start to yield 
similar results. 

The temperature dependences of the order 
parameter— and the gap in the c direction, obtained 



FIG. 9: Temperature dependences of the tunneling parameter 
Z and the scale factor S for the LuNi2B2C-Ag point contact 
in the c direction obtained with models from Refs.^^ andl25l 
The scale factor was found from fitting over the entire curve 
and over the "wings" (Scaling max) 
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FIG. 10: Temperature dependence of the order parameter 
A and energy gap Acp- and the BTK energy gap— for 
the LuNi2B2C-Ag point contact in the c direction. The T- 
dependence of the largest gap (order parameter) was esti- 
mated from fitting over the entire curve and over the "wings" . 



by the entire-curve and the "wing" BTK 2 - fittings, are 
shown in Fig.EH along with the BCS curve. As in the 
case of the scaling factor S (Fig.^J, deviation starts 
at T < 11 K, while at higher temperature both fitting 
procedures result in the same BCS-like dependence. 

The temperature dependences of the order 
parameter— and the gap; 18 i 25 for the ab plane con- 
tact are shown in Fig.^2 The calculation was made 
both over the entire curve and the "wings". It is seen 
that the curves coincide in a wider temperature interval 
and follow the BCS-dependence down to about 8 K. 
Thus, the gaps in the ab plane and c direction differ not 
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FIG. 11: Temperature dependence of the order parameter A 
and energy gap Ao (calculation by Ref. llSt) and the BTK en- 
ergy gap (calculation by Ref. l25l) for the LuNi2B2C-Ag point 
contact in the ab plane. The T-dependence of the largest gap 
(order parameter) was estimated from fitting over the entire 
curve and over the "wings". The scatter of the points at 
T — 1.5 — 4.2 K is caused by the contact instability. The 
contact resistance becomes stable above 4.5 K. 



only in magnitude but in temperature dependence as 
well. 

Let us consider in more detail the recent STM investi- 
gation of the superconducting gap anisotropy in Lu and 
Y compounds by means of STM. 19 The tunneling spectra 
obtained in this work appeared to be impossible to fit in 
terms of the traditional BCS theory. For the curve shown 
in Fig. 2 of the cited paper, the zero value of the differen- 
tial conductivity is seen up to 0.8 mV and the maximum 
of dl/dV, which corresponds approximately to the energy 
gap, is located at 2.3 mV. At the same time, the shape 
of the curve is permanent along the surface area much 
greater than the coherence length. In Fig. 2 b, at energies 
slightly less than 2 mV, there is a shoulder, which could 
be ascribed to a two-gap spectrum. Unfortunately, this 
feature is not discussed in the cited paper. The authors 
notice that the use of the model, where the broadening of 
the BCS density of states is caused by the finite life time 
of quasiparticles, implies a non-zero density of states (or 
differential conductivity) at infinitesimal bias which is at 
odd with the experimental result. Because of that, the 
authors of Ref.[l3 use a modified density of state modeled 
by the energy gap with the Gauss distribution centered 
at A with the width e. By investigation of different parts 
of the crystal surface, the authors of the cited paper ob- 
serve the correlation between the supposed anisotropy of 
the gap expressed in e/A units and the local transition 
(critical) temperature at the given spot on the crystal 
surface. With increase of the critical temperature under 
the contact, both the absolute value of the gap and its 
anisotropy increase. Such a behavior corresponds to the 
s-wave pairing, since in terms of that model the elastic 



scattering leads to a decrease of T c and isotropy of the 
gap. 

In our experiment the critical temperature coincides 
for both directions with the bulk value, as one could 
expect in a crystal with undisturbed surface. Our pa- 
rameter T allows one to estimate quantitatively the gap 
anisotropy and is similar to the parameter e/A in Refill 9l 
Namely, Y is about 1.5 times greater in the c direction 
compared with that in the ab plane. The temperature 
interval where our gap dependences deviate from BCS 
curve is also about 1.5 times larger in the c direction 
compared with that in the ab plane. Hence, we may con- 
clude that the anisotropy of the superconducting gap in 
the c direction is noticably greater than in the ab plane. 

Note that neither of the techniques^^^ can describe 
adequately the presumed situation because in both cases 
the discrepancy is determined by the distribution of the 
gap over the Fermi surface. Taking into account the poor 
efficiency of the one-gap fit at low temperatures (Fig. 5) 
and the deviation of the fitting curve at the minima of 
dV/dl for the two-gap fit (see Figs. 6, 8 which give ap- 
proximately the same discrepancy factor F as for the low- 
est temperature in Fig. 5), we can assert that the super- 
conducting gap varies continuously over the Fermi sur- 
face. It is therefore most reasonable to describe such 
curves in terms of a number of parallel-connected point 
contacts having different gaps. Their contribution to the 
total conductivity can qualitatively account for the re- 
gion of the Fermi surface containing a particular gap, as 
was considered in Ref.l27l 

Let us illustrate this approach by fitting the experi- 
mental curves for the c- and afe-directions at the low- 
est temperatures, where the discrepancy is the highest 
(Figs. 6, 8). Just as the two-gap fit JSJ, the multigap fit 
is expressed by 



dV 



dl §(A 1 ,Z)K 1 



§(A n ,Z)K n 



(6) 



where K\ + K-2 + . . . + K n = 1, Ki is proportional to 
the part of the Fermi surface with the gap Aj, and S is 
the scaling factor. The broadening parameter (7 or Y) 
was taken to be zero, and the barrier factor Z was the 
same for all the contributions. In this approximation the 
approaches of Refs.0 and|2^ are the same, and the or- 
der parameter coincides with the energy gap. As a first 
approximation, it is supposed that the coefficients Ki are 
located in the curve which is a superposition of two peaks 
(Fig. 12). The maxima of these peaks are at energies A/, 
and A p , and their slopes are described by asymmetri- 
cal Gaussian distributions with a width at half height 
of eri,<72,cr3 and 174. The heights of each peak are hi 
and hi, respectively. By varying all the fitting parame- 
ters (Aj, hi, o~k for i = 1, 2 and fc = 1 — 4) we minimized 
the average r.m.s. deviation. This procedure includes 
manually offsetting dots of the K(A) fitting curve in the 
intermediate steps of the adjustment. 

The result is shown in Figs. 13 and 14 for the c- and ab- 
directions, respectively. One can see that the multigap fit 
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TABLE I: The fitting results for superconducting parameters along two mutually perpendicular directions, using Ref.Jlg and 
the BTK theory^ Ai, A2 are the order parameters in the two-gap fitting; whereas, the subindex corresponds to the energy 
gaps in the model of Ref.llSt The pair breaking parameters 71,2 and the broadening parameter I\2 correspond to Refs.fla 
and I2EI. respectively. Z is the tunneling parameter. The temperature corresponds to that of the experiment. The relative 
contributions of each superconducting order parameter and gap are given in percent. A and F are given in meV, 7 and Z are 
dimensionless. 

Two-gap fit One-gap fit 

Ref.rS BTK2£ 
ab c ab c ab 

50% 3 62% 2 66 2 ' 74 2.25 2.55 





Ref.rS 

c 




ab 




BTK 25 
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A 1 


3.2 


40% 


3.03 


60% 


2.65 


A01 


2.9 




2.92 






A 2 


1.82 


60% 
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40% 








1.7 


A 02 


1.65 




1.82 






7i, r i 


0.016 




0.004 




0.4 


72, r 2 


0.016 




0.012 




0.4 


z 


0.558 




0.765 




0.59 


T(K) 


1.47 




1.48 




1.47 



1.88 2.34 



50% 2.2 38% 



""' 0.094 0.031 0.655 0.434 
0.14 

0.77 0.55 0.745 0.6 0.804 

1.48 1.47 1.48 1.47 1.48 



down to small values, which agrees with the {s+g) model 

of the gap nodes in some directions^ 
0.40 1 1 ° * 




Voltage (mV) 



FIG. 12: Modeling of the anisotropic superconducting energy 
gap distribution over the Fermi surface with multiple parallel 
junctions of different gaps. The contribution of a particular 
energy gap correlates with its part over the Fermi surface. 
The distribution of a particular contribution to the total con- 
ductivity is described by the Ki(Ai) dependence. The shape 
of the distribution Ki(Ai) is simulated by superposition of 
two peaks. Each peak is described by the following parame- 
ters: the energy of the central maximum A^ or A p , half- width 
of the right and left slopes 01,02 or 0-3,0-4, and the heights 
hi and /12, correspondingly. By variation of these 8 param- 
eters we attempt to reach the best coincidence between the 
experimental and theoretical curves. 



with 8 fitting parameters and a continuous distribution of 
the gaps matches the experimental curve very well. One 
should not trust this fitting literally, but it shows that the 
width of the gap distribution can be quite appreciable, 



FIG. 13: Approximation of the experimental curve by the 
use of the multigap model for c direction. T = 1.47 K, Z = 
0.559. Inset: the contribution to the total conductance from 
different parts of the Fermi surface with different gaps. The 
distribution is modeled by 31 equidistant points, A m i n =0.7 
meV, A m ax~4: meV. The resulting curve exibits two peaks 
with maxima at 1.8-2 meV and 2.5 meV. 



IV. CONCLUSION 

In this study the anisotropy of the superconducting en- 
ergy gap was measured for the first time on LuNi 2 B 2 C 
in the ab plane and c directions. It is found that at 
low temperatures the experimental curves should be de- 
scribed assuming a gap distribution over the Fermi sur- 
face. Within the two-gap model, the largest contribution 



9 




Voltage (mV) 

FIG. 14: Approximation of the experimental curve by the 
use of the multigap model for the ab plane. T = 1.48 K, 
Z — 0.755. Inset: the contribution to the total conductance 
from different parts of Fermi surface with different gaps. The 
distribution is modeled by 31 equidistant points, A m i„ = 1.0 
meV, A maa; =3.35 meV. The resulting curve exibits two peaks 
with maxima at 2 meV and 3.1 meV. 

to the total conductivity is made by the Fermi region 
with a smaller gap in the c direction and by the Fermi 
region with a greater gap in the ab plane. At T ~ 1.5 K 
the largest and the smallest order parameters in the c 
direction are A max =3.2 meV, A m j J1 = 1.82 meV; in the 
ab plane these are A max =3.03 meV, A m i„=2.0 meV. An 
attempt to fit the low temperature experimental curves 
shows that the gap is distributed starting from 0.8 meV. 
This corresponds well to the recent STM observations in 
Ref.0 discussed above. We have found the deviation 
of the temperature dependent gap from the BCS theory 
for both tested direction. This deviation is connected 
with the impossibility to describe the gap distribution in 
terms of the one-gap model, at least at low temperatures. 
The temperature range where the described deviation is 
observed for the c direction, is about 1.5 times greater 
(1.5 11. 5K) than that in the ab plane (1.5 -=-7.5 K). 
The broadening parameter T, allowing quantitative esti- 



mation of the degree of anisotropy, is also bigger in the c 
direction. 

The experimental results are described on the basis 
of the generalized BTK models and the Beloborod'ko 
theoryjM considering the electrical conductivity of ballis- 
tic S — c — N point contacts in the presence of an arbi- 
trarily penetratable potential barrier and allowing for the 
finite lifetime of Cooper pairs. For superconductors with 
a multiband electronic structure, the interband transi- 
tions of Cooper pairs may lead to their finite lifetime. 
Previously, a two-band model was suggested for nickel 
borocarbide superconductors^ The theory^ which ac- 
counts more accurately for the force of pair breaking, 
may explain the difference between the order parameter 
and the gap. 

Our next publication will concern the compound 
ErNi2B2C which exhibits a magnetic transition near 
T=6K. In this context a detailed analysis and com- 
parison of the two theoretical approachesi2*2^ applied 
to LuNi2B2C are of paramount importance for under- 
standing the results that can be obtained on magnetic 
ErNi 2 B 2 C. 
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